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We address recent momentum resolved radio frequency (RF) experiments on ultracold trapped Fermi gases of 
'"^K. We show that momentum resolved RF probes provide measurements of the centrally important fermionic 
spectral function. They also serve to remove ambiguity plaguing the interpretation of momentum integrated RF 
experiments by establishing a clear dispersion signature of pairing. We find that the temperature dependence 
of the spectral function is dramatic at unitarity, and, importantly, smooth from above to below throughout 
BCS-BEC crossover. This should be tested experimentally, given widespread predictions of first order behavior. 



Ultracold Fermi gases undergo BCS to Bose-Einstein con- 
densation (EEC) crossover via a Feshbach resonance which 
is tuned with variable magnetic field. Study of BCS-BEC 
crossover has been argued [IJ to be important for under- 
standing high temperature superconductivity. These ultracold 
Fermi gases and their extensions to optical lattices are viewed 
as quantum simulators of some of the most important prob- 
lems in condensed matter physics. However, due to the ex- 
treme small size of the Fermi gas cloud and the ultra low 
temperatures, the available experimental probes are limited by 
comparison with their condensed matter counterparts. Of par- 
ticular importance are measurements of the fermionic spec- 
tral functions, ^(k, w) for which angle resolved photoemis- 
sion spectroscopy (ARPES) has been so powerful in electronic 
systems. It is clear that counterpart experimental probes of 
A(k, w), in the ultracold gases are crucial to progress. 

Recent experiments on from the JILA group have 
now demonstrated that it is possible to measure these spec- 
tral functions using momentum resolved RF pairing gap spec- 
troscopy over a range of magnetic fields throughout the BCS- 
BEC crossover. There is a substantial advantage of using '"^K 
over the more widely studied ^Li since, for the usual Fesh- 
bach resonance around 202 G, there are no nearby competing 
resonances to introduce complications from final state inter- 
actions. In this Letter we address the theoretical basis behind 
these measurements and demonstrate that this momentum re- 
solved RF technique has, in principle, the same capabilities 
as that of ARPES. We then show that this technique can be 
used to shed light on the recent controversy [3] about the ori- 
gin of double peak structure seen in previous iQ] momentum- 
integrated RF spectra, as to whether it has to do with pairing 
or bound state effects. Both our calculations and the mea- 
surements in Ref. of the momentum resolved RF spectral 
weight support the interpretation that a generic double peak 
structure in "'"K is associated with paired atoms near the trap 
center and unpaired atoms near the trap edge. 

Importantly, we find that, near unitarity, the spectral func- 
tion has a significant temperature variation. As in ARPES 
measurements of the temperature dependent spectral function 
these experiments can provide evidence for the presence of a 
normal state gap or pseudogap above Tc at unitarity. Indeed, a 



central focus of the cuprate literature has been to address the 
evolution of angle resolved photoemission experiments from 
above to below Tc- This focus provides the motivation for the 
present work, which addresses the counterpart to these exper- 
iments in the cold Fermi gases. Our work is based on a for- 
malism that has been successfully applied to to ARPES mea- 
surements [5J in the high temperature superconductors. The 
corresponding physical picture of the pseudogap which opens 
at a temperature T* > Tc is that it arises from a stronger-than- 
BCS attraction leading to preformed, normal state pairs. The 
famous BCS paper highlights an additional challenge which 
must be met by theory: namely, explaining a smooth or sec- 
ond order phase transition at Tc- In the presence of BCS-BEC 
crossover we find that the transition should also be second or- 
der. This is at odds with widespread predictions of first 
order behavior in the literature. 

In addressing RF spectroscopy we can ignore final state ef- 
fects i^l for the particular resonance under consideration. In 
the recent JILA RF experiments |2] the momentum of state 3 
atoms is measured using time-of-flight imaging, in conjunc- 
tion with 3D distribution reconstruction. Since the 3D gas in 
a single trap is isotropic, detailed angular information is irrel- 
evant. In a homogeneous system, the contributions to the RF 
current from state 2 to state 3 are of two types corresponding 
to transitions associated with the breaking of 1-2 pairs in the 
condensate or with the promotion of already existing thermal 
fermionic excitations to level 3. In RF transitions, a domi- 
nantly large fraction {Q. l) of the photon energy hv is used to 
excite the atoms in hyperfine level 2 to another internal state 
of the atom, coiTesponding to level 3, where 51^ is the photon 
energy needed when the atoms in level 2 are free. By contrast, 
in the ARPES case, there is an "emission" process in that the 
photon energy is converted into the electron kinetic energy so 
that the electrons are knocked out of the crystal. 

In the many-body theoretical derivation, it is convenient to 
choose the Fermi level of state 2 as zero energy, so that we 
define il^ as the overall energy level splitting between the 
chemical potentials for atoms in level 3 and level 2. When 
pairing is present, the Green's function G{K) for this state is 
gapped, and the Green's function Gg"' [K) for state 3 is free. 
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Figure 1: Contour plot of the occupied spectral intensity 
A{k,Lu)f{Lo)k'^/2-ir'^ at unitarity in the homogeneous case at 
T/Tc ~ 1.9. The three lines in the contour plots correspond to 
the two quasiparticle dispersions iiJ]^ and the free atom dispersion, 
respectively, where we take fi = 0.7Ef, A = 0.5Ef, T = OATf, 
and 7 = O.I-Bf- 



We then have for the momentum-resolved RF current 

\Tk\' 



I{k,diy) = J dLu'A{k,u}' -n)A3{k,uj') 



(2^)2 

\Tk\ 



2-K 



(1) 



where K = (k, iiOn) is the four-momentum, and Tk is the 
transition matrix element, which we will take to be a constant 
and then set to unity. This expression has been previously 
derived in Ref. [9], except there was also a momentum inte- 
gral. We have also assumed that state 3 is essentially empty 
[/(w' + ri'^) — 0], and has a free dispersion ^3(k, lj) = 
27r(5(ti; — ek+A*3~M)- Note here = fi— /x+/i3 = hv—flL is 
the RF detuning and uj = e]^ — Siy is the energy of state 2 mea- 
sured with respect to the Fermi level, where Ck — /2m — fi, 
jjL and are the chemical potentials of atoms in state 2 and 3, 
respectively, and we take h—\. The last line of Eq. ([TJ is for- 
mally identical to the expression for ARPES signal intensity. 
In this way, it should be clear that momentum-resolved RF 
is equivalent to ARPES. Both measure the fermionic spectral 
function. 

Beginning with Refs. [l^ and n, there have been studies 



of the spectral function associated with BCS-BEC crossover 
which show that, beyond the BCS regime, there is a spec- 
tral gap or pseudogap. Related work has also been performed 
by other groups |12]. Detailed numerical calculations 1,1 1,] 
of the related normal state self energy (called Spg) have led 
to a rather simple BCS-like form, which has been shown 
phenomenologically Iil3i] to be quite successful for describ- 
ing the normal phase of the cuprates. Following our earlier 
body of work [Sl S [13] we presume that the self energy for a 
short-range interaction consists of a term arising from the con- 



densed, Yisc, as well as noncondensed pairs (E^g). We write 
i;(k, uj) = Spg(k, bj) + EscCk, Lo), so that 



Epg(k,w) 
Ssc(k,a;) 



A2 

pg 



W + Ek 



(2) 
(3) 



Here 7 7^ reflects a phenomenological broadening term, as- 
sociated with a finite life time for incoherent pairs, and Eo 
represents an "incoherent" background contribution, includ- 
ing that from the particle-hole channel. Ssc is associated with 
long-lived, condensed Cooper pairs and so it is of the same 
form as Epg but without the broadening. With this self en- 
ergy in the Green's functions, the resulting spectral function, 
^(k, = — 2ImG(k, + iO) can be readily determined. 
There is, moreover, a well defined procedure for computing 
the two contributions (from the condensed and noncondensed 
pairs) to the excitation gap which we do not repeat here. This 
formalism has been applied in the cold Fermi gas ex peri ments 
lilllslEIlT'] as well as in the cuprate Hterature lUllsll. 

In Fig. 1, we present an intensity map as a function of 
the single particle energy and the wave vector k for the states 
which contribute to the RF current in a homogeneous system. 
The yellow, red regions of the figure indicates those pairs of 
points {E, k) where the greatest number of contributing atoms 
lie. The temperature here is chosen to be relatively high, 
around l.OTc in order to have appreciable contributions to the 
RF current from pre-existing thermal fermionic excitations. 
The intensity map indicates upward and downward dispersing 
contributions. These correspond, to a good approximation, to 
the two RF transitions to state 3 from state 2 with dispersions 
(Ey^ + fi) and {—E^ + ji) relative to the bottom of the band. 
Here E^ ~ y^eJT+A^ is the quasiparticle excitation energy, 
and = A^^ + A^^ defines the total gap A, with Asc = 
above Tc- The width of the peaks in this contour plot comes 
exclusively from the incoherent terms 7 and Sq. One can see 
that the bulk of the cuiTent even at this high temperature is as- 
sociated with the pair states which are broken in the process of 
the RF excitation. This figure describes in a conceptual way, 
how this intensity map can be used to compare with a broad- 
ened BCS-like form for the spectral function, which fits very 
well the two branches shown in the figure corresponding to 
upward and downward dispersing curves. From the intensity 
map, one can determine the spectral gap by fitting the spectral 
function or the dispersion to a broadened BCS-like form, as 
has been done in Ref. |2]. 

We stress that the above calculations are for the homo- 
geneous case and it is important to extend them to include 
the effects of a trap, where the density and gap reach max- 
imum at the trap center and vanish at the trap edge. This 
is done within a local density approximation (LDA). Once 
the trap is included the simple analogy between the elec- 
tronic ARPES experiments and momentum resolved photoe- 
mission spectroscopy is invalidated. While 6v is constant 
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Figure 2: Contour plots of the occupied spectral intensity A(k,Lu)f{Lo)k'^/2iT^N at unitarity for (a) T/Tc = 1.35, (b) 0.8, and (c) 0.1, from 
left to right. Here Eo = 0.25E% and 7 = 0.25{T/Tc)Ep at the trap center, and both scale with local density as r varies, Tc ~ 0.27Tp. 
The results were convoluted with a Gaussian broadening curve of width a — Q.22E%, whose value is taken from Ref. |2]. As T decreases, 
the spectral weight shifts from a free atom peak (upward dispersing curve) to a paired atom peak (downward dispersing curve). Here the lines 
indicate the loci of the peaks in the EDCs. 



across the sample, the single particle excitation energy uj — 
k'^ /2m — ^{r) — Si/ will necessarily bear the inhomogeneity 
effect of the spatial dependence of /i(r); after time of flight 
imaging, one cannot distinguish atoms from different radii r. 
It is thus more convenient to measure the single particle en- 
ergy from the band bottom, uj + ^{r) ~ k'^ /2m — 5v, which 
is constant across the trap for given k. Nevertheless, many of 
the central features survive. While in a homogeneous system 
there will be two branches as described above, in a trap, there 
is a third branch which appears as well. This corresponds to 
essentially free atoms at the trap edge which will contribute 
significantly to the RF current. This contribution dominates 
that of the thermally excited quasiparticles at both low and 
high T. It is this branch which is also upward dispersing 
and makes it rather difficult to see the contributions from pre- 
existing thermally broken pairs. With improved experimental 
resolution, this momentum resolved RF probe can in principle 
distinguish these two contributions by their dispersion. 

The results of the LDA-based calculations for the intensity 
map plots are shown in Fig. 2 for the unitary gas in a trap, for 
different temperatures, calculated with Sq = 7 = Q.2bE^, 
and convoluted with an experimental Gaussian line of width 
cr = Q.22Tp. The vertical axis represents the trap averaged 
single particle energy, {oj + fj,)/Ep, where Ep = {kp)'^/2m 
is the global Fermi energy in the non-interacting limit. The 
temperatures are T/Tc = 1.35 (above Tc), 0.8 (slightly be- 
low Tc), and 0.1 (well below Tc), from left to right. At the 
highest T the central notable feature is a single upward dis- 
persing curve which fits the free particle dispersion; the spec- 
tral weight is dominated by contributions from atoms at the 
trap edge. This dispersion can be readily differentiated from 
that associated with pre-existing thermally broken pairs which 
varies as E^^ + n and of course, depends on the distribution of 
energy gaps A(r). 

As the temperature decreases a second (downward dispers- 
ing) branch becomes evident. The spectral weight shifts to the 
lower branch gradually. The intensity map first bifurcates and 
eventually becomes dominated by the lower branch at very 



low T, where essentially all atoms are paired. Here the lower 
branch is associated with the breaking of a condensate pair. 
It, moreover, contains trap averaging effects. In the absence 
of final state effects, this can be used to determine the pair- 
ing gap. The separation of the two peaks can be difficult to 
discern until sufficiently high values of k. 

That there are no abrupt changes at the superfluid transi- 
tion distinguishes this theoretical scheme from others (e.g., 
Refs. iaLJ) in the literature. Just as for the high temperature 
superconductors, there are no first order transitions as the sys- 
tem evolves from above to below Tc. The behavior at unitarity 
is in some sense smoother than in the strict BCS theory, since 
we find that a pairing gap is present above Tc, which shows 
up as a bifurcation in the intensity plots. There are indications 
120 from the experimental momentum resolved RF plots that 
the behavior evolves smoothly from the normal to the super- 
fluid phase, but detailed quantitative analyses to rule in or out 
possible first order transitions are lacking. 

In Fig. 3A we plot the energy distribution curves (EDC) 
for a series of momenta fc in a unitary trapped Fermi gas at 
T/Tc = 1.1 as a function of single particle energy. Here 
we have taken a somewhat larger intrinsic and instrumental 
broadening than in Fig. 1. These parameters are seen to op- 
timize semi-quantitative agreement with the data. Thus the 
EDC in Fig. 3A can be compared with their counterparts plot- 
ted in Fig. 4 in Ref. |2]. At small k there is a single peak 
which separates into two at higher k values. The data are not 
sufficiently smooth to be assured that two peaks are visible at 
larger fc, but they are consistent with this picture. As before, 
we attribute the separation at higher k to the fact that the up- 
ward and downward dispersing curves become most well sep- 
arated in this regime. Better agreement is expected with im- 
proved signal-to-noise ratio in the data. Moreover, including 
detailed incoherent background contributions in the fermionic 
self-energy may also improve the agreement with experiment. 

Figure 3B presents the corresponding spectral intensity 
map. The white dashed curve represents a fit of the experi- 
mental peak dispersion [2] while the solid white curve is the 
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Figure 3: (A) Energy distribution curves for a series of momenta 
A: as a function of single particle energy, and (B) corresponding 
occupied spectral intensity map, in a unitary trapped Fermi gas at 
T/Tc = 1.1. Here Eo = 0.35E°p and 7 = 0.3SE°p at the trap cen- 
ter, and a — 0.3Ep. In (B), the red curves represent the free atom 
dispersion, while the white solid and dashed curves are the quasi- 
particle dispersion obtained theoretically and experimentally |2], re- 
spectively, via fitting the EDCs with a single Gaussian distribution. 

theoretical counterpart. Here, as in experiment we have fit the 
EDC to a single Gaussian peak. The comparison between the 
two white curves shows semi-quantitative consistency. More- 
over, both the solid and dashed white curves can be well fit to 
the BCS dispersion involving E^^, as was originally proposed 
in Ref. |10]. Indeed Fig. 3B seems to capture the essential 
results shown in Fig. 3a in Ref. [2]. With higher experimental 
resolution it should be possible to obtain more direct informa- 
tion about the mean gap size. In these experiments on ^"K, 
final state effects should be of less importance in complicat- 
ing the interpretation. The fact that the experiments were done 
near Tc suggests that there is a sizable pseudogap at and above 
Tc at unitarity. 

We turn finally to the BEC side of resonance in Fig. 4. Here 
the intensity map is presented for l/k^a = 1.1 at a tempera- 
ture T/Tc = 1.3, where w 0.3T^. Even at this relatively 
high temperature there is virtually no weight in the free atom 
peak, associated with atoms at the trap edge. All the atoms 
are bound into pairs and there is only one downward dispers- 
ing peak in the plot. This behavior can be contrasted with that 
shown in Fig. 3c of Ref. |2] where there is a free atom con- 
tribution, which has been attributed to unpaired atoms out of 
chemical equilibrium with the pairs. In addition the down- 
ward dispersing peak in the left panel is rather sharp. In con- 
trast, the large peak width in Ref. [2] seems to suggest inco- 
herent contributions to the fermionic self energy, which may 
have to do with chemical non-equilibrium. The underlying 
theory here presumes that there are still well defined quasi- 
particles on the BEC side of resonance. We present in Fig. 3b 
the case in which we have considered a substantially larger in- 
trinsic broadening 7 « A to represent some of the qualitative 
features of the experimental data. 

In conclusion, we have shown that momentum resolved RF 
spectroscopy can be used to measure the centrally important 
spectral function of ultracold Fermi gases. Other spectro- 
scopic experimental tools have been proposed in the literature 
fl9il . although currently there are no counterpart experimental 



Figure 4: Occupied spectral intensity maps at l/k%a = 1.1 in the 
BEC case and T /Tc = 1.3. The broadening parameters at the 
trap center are Eo = O.SSiSj- and 7 = QA^E% for the left panel 
and Eo = l.OEp and 7 = l.iE^p for the right panel. For both, 
a = 0.25iJ^. Coherent sharp quasiparticle peaks can be seen in the 
left panel. Larger broadening in the right panel makes the coherent 
peak less pronounced, mimicking possible effects of chemical non- 
equilibrium. The red curve represents free atom dispersion, and the 
solid and dashed white curves represent the theoretical and experi- 
mental quasiparticle dispersion, respectively. 

studies. We show that our calculated spectral intensity maps 
are in semi-quantitative agreement with experiments on ^°K. 
Because they establish the actual dispersion associated with 
the RF current contributions, these experiments remove am- 
biguity associated with whether pairing contributions might 
be confused with final bound state effects fi] . Both our the- 
ory and the experiments in Ref. reveal that a double peak 
structure reminiscent of that seen in previous RF spectra \^ is 
indeed associated with paired atoms at the trap center and un- 
paired atoms at the trap edge. At unitarity, we show that there 
is substantial temperature variation in the spectral weight dis- 
tribution. Theoretically, the spectral intensity maps and quasi- 
particle dispersion evolve smoothly with temperature across 
Tc- Whether first order transitions exist, as some have pre- 
dicted, will require further experiments. 
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